Sub-Alfvenic Non-Ideal MHD Turbulence Simulations with Ambipolar Diffusion: I. Turbulence Statistics

Most numerical investigations on the role of magnetic fields in turbulent molecular clouds (MCs) are based on ideal magneto-hydrodynamics (MHD). However, MCs are weakly ionized, so that the time scale required for the magnetic field to diffuse through the neutral component of the plasma by ambipolar diffusion (AD) can be comparable to the dynamical time scale. We have performed a series of 256^3 and 512^3 simulations on supersonic but sub-Alfvenic turbulent systems with AD using the Heavy-Ion Approximation developed in Li, McKee, & Klein (2006). Our calculations are based on the assumption that the number of ions is conserved, but we show that these results approximately apply to the case of time-dependent ionization in molecular clouds as well. Convergence studies allow us to determine the optimal value of the ionization mass fraction when using the heavy-ion approximation for low Mach number, sub-Alfvenic turbulent systems. We find that ambipolar diffusion steepens the velocity and magnetic power spectra compared to the ideal MHD case. Changes in the density PDF, total magnetic energy, and ionization fraction are determined as a function of the AD Reynolds number. The power spectra for the neutral gas properties of a strongly magnetized medium with a low AD Reynolds number are similar to those for a weakly magnetized medium; in particular, the power spectrum of the neutral velocity is close to that for Burgers turbulence.Comment: 37 pages, 11 figures, 4 table

The Formation and Evolution of Wide-orbit Stellar Multiples In Magnetized Clouds

Stars rarely form in isolation. Nearly half of the stars in the Milky Way have a companion, and this fraction increases in star-forming regions. However, why some dense cores and filaments form bound pairs while others form single stars remains unclear. We present a set of three-dimensional, gravo-magnetohydrodynamic simulations of turbulent star-forming clouds, aimed at understanding the formation and evolution of multiple-star systems formed through large-scale (greater than or similar to 10(3) au) turbulent fragmentation. We investigate three global magnetic field strengths, with global mass-to-flux ratios of mu(phi) = 2, 8, and 32. The initial separations of protostars in multiples depend on the global magnetic field strength, with stronger magnetic fields (e.g., mu(phi)= 2) suppressing fragmentation on smaller scales. The overall multiplicity fraction (MF) is between 0.4 and 0.6 for our strong and intermediate magnetic field strengths, which is in agreement with observations. The weak field case has a lower fraction. The MF is relatively constant throughout the simulations, even though stellar densities increase as collapse continues. While the MF rarely exceeds 60% in all three simulations, over 80% of all protostars are part of a binary system at some point. We additionally find that the distribution of binary spin misalignment angles is consistent with a randomized distribution. In all three simulations, several binaries originate with wide separations and dynamically evolve to less than or similar to 10(2) au separations. We show that a simple model of mass accretion and dynamical friction with the gas can explain this orbital evolution.This item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at [email protected]